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Abstract 

It is pointed out that the internal structure of the narrow resonance X(3872) at the 
D^D*^ threshold can be studied in some detail by measuring the rate and the spectra 
in the decays X(3872) D^D^tt^ and X(3872) D^D^-i. In particular, if this 
resonance contains a dominant 'molecular' component DD* it DD*, this component 
can be revealed and studied by a distinct pattern of interference between the underlying 
decays of D*^ and D*^ whose coherence is ensured by fixed (but yet unknown) C parity 
of the X(3872). 



The recently observedj^^ by the Belle Collaboration narrow resonance X(3872) at 3872. 0± 
0.6±0.5Me1/ (and confirmed by CDF[21 at 3871.4±0.7±0.4 Mel/) decaying into vr+vr-J/?/' 
is within 0.2 ± 0.7 MeV from the D^D*^ threshold. The extreme proximity of the resonance 
to the threshold naturally invites the suggestion [21 E] that its wave function may have a 
large component with a pair of neutral pseudoscalar and vector D (D) (anti)mesons. The 
spatial separation of the mesons in this component is sufficiently large for the mesons to 
retain their individual structure. Such configuration would clearly realize the long-standing 
conjecture 13 EI of existence of "molecular charmonium" i.e. of resonances, which essentially 
are loosely bound states of charmed hadrons. 

Clearly, the approximately 1 MeV or less scale for the energy gap w between the resonance 
and the threshold is quite likely to result in a completely different weight in the wave function 
of the X(3872) of the pairs of neutral and charged D mesons, since the threshold for the 
charged D^D*^ pairs is another 8.1 ± 0.5 MeV higher, which is 'far' in the scale of w. Thus 
the isospin is likely to be strongly violated in the X(3872) resonance, which in the case 
if X(3872) is even under C parity, would allow the observed decay X(3872) vr+vr" J/?/^ 
to be in fact occurring due to the decay X(3872) — > p^Jip E], in agreement with the 
very strong peaking at the maximal value of the spectrum of the invariant mass of the two 
pionsj^. One can trivially notice that this conjecture can be readily tested by a search for 
decay involving neutral pions: X(3872) —* it^Ti^J/ip. If indeed the pions emerge from the p° 
resonance, the process with neutral pions should be absent. The decay with neutral pions 
is also forbidden, and the observed decay X(3872) — > -K^-K^J/ip is manifestly due to an 
1 = 1 component of X(3872), in the general case of C(X) = +1, even if the p° dominance 
is not confirmed. Alternatively, if the discussed resonance is a C = —1 state, the decay 
X(3872) -K^-K^J/ip is allowed with the dipion being in the / = isospin state, and the 
relation r(X tt^tt~J/iIj) = 2r(X — > Tr^TT^J/ip) should hold to a good accuracy. (Any 
significant presence of an / = 2 state of the dipion would obviously be totally exotic.) 

Naturally, further study of the properties of X(3872) will likely involve other possible 
decays of this resonance, including the decays related to the underlying transitions D*^ —>■ 
D^TT^, D*^ -D°7, and the corresponding transitions between the anti- mesons |Sl !!]• The 
main purpose of the present paper is to point out that at the characteristic momenta of the 
mesons in the wave function of the molecular D^D*^ component of the X(3872) resonance the 
parameters of these decays should likely be measurably different from those of an incoherent 
sum of decays of free D*^ and Z)*° mesons. Rather the rates and the spectra of the decays 
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X(3872) D^D^TT^ and X(3872) D^D^^i should exhibit binding effects and a significant 
interference between the underlying decays of vector mesons and anti-mesons. Thus an 
experimental study of these decays may reveal rather fine details of the structure of the 
X(3872) resonance. In other words, the Dalitz plots of these decays would provide a "CAT 
scan" of the actual wave function of the mesons inside X(3872). 

It can be noted, that the spatial (momentum) dependence of the main part of the wave 
function of the mesons can be described, in a way, similar to that used for deuteronjSj, and the 
essential unknown parameter for this part is the overall normalization, which represents the 
weight of the molecular component in the wave function of the X(3872) resonance. Indeed, 
at the gap energy w"^! MeV the dynamics of the D^D*^ [D^D*^) meson pair is determined 



by momenta of order k = y2/i|e| ^4:5 MeV, where /i ~ 966 MeV is the reduced mass of 
the system made of pseudoscalar and vector neutral D mesons. Thus the characteristic 
distances are far beyond the range of the strong interaction, and the wave function at 
those characteristic distances is in fact given by the Schrodinger equation for free motion. 
On the other hand the value of k may well be comparable with the momentum p of the pion 
emitted in the decay X(3872) D^D'^tt^ {pq = 43 MeV for a decay of a free D*^ meson), 
which would give rise to large binding and interference effects in the decay. In the case of 
radiative decay X(3872) — > D^D^'j, the representative value of the photon momentum k is 
that in a free D*^ decay: ko = 137 MeV. Although this value looks large as compared to k, 
it will be shown that the relative magnitude of the interference effect is determined by the 
expression {2k,/ k) arctan(/c/2K) and is significant for this decay as well. 

In the following discussion it is assumed for definiteness that X(3872) is below the D^D*^ 
threshold, so that w = m + M — M{X) is a positive quantity, where m = M{D^) ^ 
1864.5 MeV and M = M{D*^) ^ 2006.7 MeF. A generalization to the case where X(3872) 
is just above the threshold can be done by analytical continuation. Also for definiteness it is 
assumed here that the mesons inside the X(3872) are in the S wave, which quite plausibly is 
the actual situation. This obviously corresponds to J^'"{X) equal to either l^"*", or 1^ . If 
it further turns out that the quantum numbers of X(3872) are different, the orbital motion 
of the mesons can readily be accounted for by a straightforward modification of the formulas 
presented below. Thus the wave function of the relative motion of the mesons is considered 
here as given by the standard S wave expression 





(1) 
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where is the overall weight of the considered here molecular component in the X(3872) 
resonance. For a purely molecular system = 1, while realistically one would expect < 1 
thus allowing for some admixture in the wave function of X(3872) of other states (e.g. cc, 
D+D*-, etc.). 

The amplitude for the decay D*^ D^tt^ can be written as 

Ad-dtt = gi^-p) , (2) 

where e is the polarization amplitude of the D* and p is the momentum of the pion. The 
coupling constant g is related to the width = r{D*^ —>■ D^tt^) as^ F^r = Ifi'^IPo/Svr. The 
rate F^^ can be estimated from the isotopic symmetry and the known[7j total width of D*^ 
(96 ± 22KeV) and the branching ratio B{D*+ D+tt^) = (30.7 ± 0.5)%, and also taking 
into account the slight difference in the kinematics: F^ = 43 ± lOKeV. In terms of the 
coupling g this leads to a quite reasonable estimate \g^^\ = 315 ± 36MeV. 

For a system of a vector and a pseudoscalar mesons, with a definite C parity rj the 
amplitude of decay into D^D "vrO is contributed by both the decay D*^ — > D^tc^ and its 
charge-conjugate D*^ Z)°7r°. Taking into account that (^(vr*^) = +1, and performing the 
standard transition to the center-of-mass coordinate R and the relative coordinate r the 
amplitude of the decay X(3872) D^D^n^ can be written as 

{D\q^)D\q,)n'{p)\Hn^nn\X{e,P = 0)) = 

(27rf5^'\qi + q2+p)gie-p) [^(gl) + r/ 0(gi)] , (3) 

where qi (^2) is the momentum of the final D (D) meson, p is the pion momentum, and P 
is the momentum of the initial X(3872) resonance, which is set to zero, corresponding to 
consideration in the rest frame of the X(3872). (In eq.Q a use is made of the momentum 
conservation relation in this specific frame: qi + q2 + p = 0, which somewhat simplifies the 
formula and the subsequent discussion.) Finally, (f){q) is the wave function of the relative 
motion in the momentum representation: 



0(g) = J il,*{f)e"^-U\ . (4) 
The 'free' wave function in eq.(^ in the momentum space reads as 



g"^ + 

^The nonrelativistic normalization is used here for the wave functions of the D mesons, but not for the 
pion. 
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It should be noted, that in the expression in eq.Q it is assumed that the final D and D 
mesons move as free particles, i.e. that the wave function of each is a plane wave exp(zg ■ r). 
Such assumption looks quite reasonable, since the final D mesons are produced at large 
distances of order from each other, i.e. beyond the range of strong interaction. This 
behavior would be invalid if there were a resonance or a bound state of the pseudoscalar 
D mesons very close to their threshold, similar to the X(3872) state at the threshold of 
D^D*^. Existence of such resonance would certainly be a new phenomenon by itself, and 
would require a separate consideration. Here it is assumed that no singularity exists in the 
spectrum of DD pairs within at least few MeV near their threshold. 

Using eq.Q the expression for the decay rate can be written in terms of (f){q) in the 
textbook form: 

dT{X ^ D°D°7r°) = 

2 

= \9^\ \<l>{q2) + V <I>{Qi)\^ dql dq^ , (6) 
where A = M — m = 142.12 ± 0.07 MeV is the difference between the masses of the vector 



and pseudoscalar neutral D mesons, and E^^ = yp^ + is the energy of the pion. The 
intermediate expression with un-integrated delta-functions is convenient for discussing the 
limiting case of loose binding, k 0, while the final one is the standard Dalitz type and 
thus is convenient for discussing the decay parameters in terms of the Dalitz plot. (Clearly, 
in the latter expression the value of is uniquely determined through the conservation laws 
by the values of qi and q2.) 

In the limit of no binding (k 0) the momentum space wave function can be replaced 
as ^'^ {2tc)^ 5^^\q) , and the intermediate expression splits into two noninterfering 

terms, corresponding to independent 'free' decays D*^ D^ti^ (g2 = 0), and Z)*° D^n^ 
[qi = 0), thus recovering the naive expression for the total width F = 2^^ F^r (and the trivial 
kinematics). However at a value of k, comparable with pq the spread of the momentum space 
wave function and the interference effects are essential. 

In order to present the results of calculation of from eq.®, we write the total rate of the 
discussed decay in the form 

F(X ^ D'D\'') = 2e F^ [Aiw) + r]Biw)] , (7) 
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where A{w) describes the incoherent contribution of the decays of individual D*^ and Z)*°, 
and B{w) describes the effect of the interference between these two processes. The resuh of 
a numerical calculation of the terms A and B with the wave function from eq.© is shown 
in Fig.l. It is seen from the plot, that the discussed effects reach quite sizeable magnitude 
starting already from small values of the binding energy w ~ OAMeV. In particular the 
interference between the two wave functions in eq.® significantly enhances the discussed 
decay if the C parity of X(3872) is positive {rj = +1) and suppresses the rate in the case of 
negative C parity {rj = —1). 




Figure 1: The non-coherent contribution A{w) (solid line) and the interference term B{w) 
(dashed) as defined in eq.Q, calculated by a numerical integration in eq.®- 



The sign of the interference term is reversed in the radiative decay X(3872) — D^D^'y 
due to the negative C parity of the photon. The general expression for the decay rate, 
analogously to eq.®, has the form 

dT{X D^D^'-i) = (8) 

I (4 \m) - .m)? ^ (a — * - 1; - 1;) (f. + * + i'^^ '^'^ t . 

where k is the momentum of the photon, is the width of the 'free' decay D*^ D^'j (from 
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the available data one can estimate ^ 26 ± 7 KeV) , and ko ^ 137 MeV is the photon 
energy in the 'free' decay. Since the binding energy w is in any case very small in comparison 
with A, one can neglect the small shift in the energy k of the photon in the decay of X(3872) 
in comparison with /cq- Furthermore, the effect of the recoil of the heavy mesons, when their 
momentum changes on the scale of both k and k^, contributes very little to the energy 
balance, and one can perform the integration over one of the heavy meson momenta by 
neglecting the kinematical constraint on it. Also making use of the normalization condition: 
/ \(j){q)\'^d^q/{2'ir)^ = one readily arrives at the following expression for the total rate 



Clearly, the interference term described by this formula is quite substantial even at very 
moderate values of n/ko: e.g. it amounts to 0.32 already at w = 0.1 MeV, and to 0.71 at 
w = 0.5 MeV. 

One might argue that the integral for the interference term in eq.Q is mainly contributed 
by the behavior of the wave function at momenta of order g ~ A;o/2 ^ 70MeV, i.e. larger 
than K. However the corresponding distances r ~ 2/ko are still beyond the range of strong 
forces, and the 'free' approximation for the wave function should still be applicable. 

As is already discussed, the free motion wave function (Q) is justified only at distances 
beyond the range of strong interactions, and thus it fails to properly describe the dynamics 
at shorter distances, r ^m~^. At those distances, i.e. in the 'core' of the system, the mesons 
strongly overlap, and the whole 'molecular' picture of individual heavy mesons is likely to 
be inapplicable. It is not known at present, how significant the non-molecular core part 
of the wave function is, and in particular what is its contribution to the amplitudes of the 
discussed decays. It is clear however that the possible core contribution to these decays 
should lead to significantly larger than k values of the momentum transfer to the final heavy 
mesons in the decays. Thus this contribution can be revealed by studying the momentum 
distribution of the D and D mesons produced in the decays. The core contribution should 
rather uniformly populate the Dalitz plot, including the events, where both heavy mesons 
recoil with a momentum significantly larger than n, up to the kinematical limits of the 




(9) 



Using the expression (0) for the momentum space wave function, one finally finds 




(10) 
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Dalitz plot. On the contrary, the discussed here 'molecular' contribution mainly populates 
the regions, where one of the heavy mesons (the spectator) has a recoil momentum of order 
K,. In other words, a study of the Dalitz plot of the discussed decays would allow to literally 
scan the internal structure of the X(3872) resonance and, possibly, to see both the molecular 
and core components of its internal dynamics. 

When this work was finished, there appeared the paper |8j, where possible properties 
of the X(3872) resonance are discussed in connection with a further study of the decay 
X(3872) 'n-^7T~ J/ip, including the possibility of this resonance being dominantly a molec- 
ular type state. 

This work is supported in part by the DOE grant DE-FG02-94ER40823. 
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